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ABSTRACT 


FI Nino is conventionally defined as an anomalous and persistent warming ol the wateis oil the 
coasts of Ecuador and Peru in the eastern equatorial Pacific, having onset usually in southern hennsphenc 
summer/fall. Examined here are some of the statistical aspects ol El Nino occurrences, especially as they 
relate to the normal distribution and to possible associations with volcanic solar, and geomagnetic 
activity For example, the elapsed time between successive onsets ol near moderate-to-stiongei 1 1 
(over the last century and a half) is about 3.7 ± 0.6 years (at the 95 percent level o confidence) having 
a standard deviation of 1.8 +0.6/-0.3 years, a range ol 1 to 8 years, and a med.an/mode va ue o . 
years The distribution of elapsed times between successive onsets appears normally distributed, 
although a positive skew is clearly seen and greater variation has occurred within the last 50 to 75 years. 

For the interval 1848 to 1987, El Nino events of moderate strength outnumber those ol strong intensi y y 

24 to 14, and, as yet, there has never been a single occurrence of back-to-back strong events. 

A comparison of El Nino onsets to major volcanic activity (where a major eruption is taken to be 
one whose Volcanic Explosivity Index VEI ^4 and/or Dust Veil Index DVI >250) results in the lollop- 
ing- About 40 and 60 percent, respectively, of the El Nino events are preceded within one or two years by- 
major volcanic eruptions in the tropics, while about 70 and 80 percent, respectively, are preceded within 
one or two years by major eruptions somewhere in the world (i.e., ignoring latitude). Also, a number ol 
El Nino events cannot be directly linked with any known major volcanic eruption that preceded them 
within three years (some within several), the most recent example being the 197--73 strong El Nino. 

A comparison of El Nino onsets to specific phases of the solar/geomagnetic cycles results in the 
following: Nearly two-thirds of the El Nino events had their onsets when annual sunspot number was 
below 54 3 (the average annual sunspot number for the interval 1848 to 1987), the bulk (about 70 pei- 
cent) of which were moderate in strength. For Cycles 10 to 22, every sunspot cycle either had an onset ol 
El Nino during the window ± 1 year from solar minimum or already had one in progress at that time 
(Cycle 13). Also, about 70 percent of the El Nino events had their onsets during the declining portion o 
the sunspot cycle (which corresponds to the interval when geomagnetic activity is usually greatest), win 
p ou t of 14 (about 85 percent) strong El Nino events having onsets during the declining portion ol the 
sunspot cycle. Lastly, evidence was found (at >95 percent level of confidence) favoring an inverse 
correlation between the number of El Nino events per sunspot cycle and the size ot the eye e which d 
true suggests that fewer than three El Nino events may occur during the present sunspot cycle (Cycle 
1986 uTca. 1996) and that more El Nino events may have occurred, on average, during the Maunder 

minimum (1645 to 1715). 

With regard to the “very strong" El Nino of 1 982-83, it is noted that, although it may very well be 
related to the 1982 eruptions of El Chichon, the event occurred essentially “on time (with respect to the 
past behavior of elapsed times between successive El Nino events; a moderate-to-stronger b Nino was 
expected during the interval 1978 to 1982, assuming that El Nino occurrences are normally distributed 
having a mean elapsed time between successive onsets of 4 years and a standard deviation ot 2 years and 
a last known occurrence in 1976). Also, although not widely recognized, the whole ot 198_ was a rect rd 
year for geomagnetic activity (based on the aa geomagnetic index, with the aa index registering an a - 
time high in February 1982), perhaps, important for determining a possible “trigger lor this and other 

El Nino events. 
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A major feature of this study is an extensive bibliography (325 entries) on El Nino and volcanic- 
solar-geomagnetic effects on climate. Also, included is a tabular listing of the 94 major volcanic 
eruptions of 1835 to 1986. 
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NASA TECHNICAL PAPER 


ON THE STATISTICS OF EL NINO OCCURRENCES AND THE RELATIONSHIP OF 
EL NINO TO VOLCANIC AND SOLAR/GEOMAGNETIC ACTIVITY 


I. INTRODUCTION 


The term "El Nino" was originally used to deseribe the annual warming of the waters off the 
western coast of South America (in particular. Ecuador and Peru) near to just lollowing Christmas. As 
such. El Nino was regarded as a purely regional and short-lived oceanic phenomenon, of little interest to 
people living outside the immediate area and ol little consequence in relation to the world s weather and 
climate. 

Within the past 25 years, however, this perception has radically changed. For example, today the 
term “El Nino” is applied exclusively to only those anomalous and pervasive warmings ot the coastal 
waters off Ecuador and Peru that are of a more persistent and intense nature, lasting typically several 
months to more than a year and occurring once every few to several years [52, 221 1. Further, changes in 
sea surface temperature (SST) in the eastern/central equatorial Pacilic, which is the primary ocean- 
ographic signal for El Nino and its counterpart, "La Nina (also called “Warm Event and Cold Event, 
respectively), have been found to often correlate with changes in the ditlerence o! atmospheric sea level 
pressure (SLP) between widely separated stations in the Pacilic Basin (typically, Tahiti or Easter Island 
and Darwin. Australia; this difference is called the "Southern Oscillation Index" or SOI and serves as the 
primary atmospheric signal lor the Southern Oscillation, a somewhat orderly rearrangement or progres- 
sion of meteorological phenomena in the region of the Indian and Pacific Oceans); hence, El Nino is now 
perceived to be part of a large-scale, dynamic, atmospheric-oceanic, climatic feedback system ol global 
proportion [1, 7-8, 12-15, 18, 20, 34-36, 39-41, 63-64, 69, 74, 86. 89. 91, 101, 106, 1 16-117, 128, 
135-136, 138. 148, 162, 173, 178. 181, 189, 192. 195, 201, 206-210, 220. 225-229, 231-234, 238. 
244-246, 249, 273-274, 286-287. 309-311, 313-316, 318-32()|. 

While it is generally accepted that LI Nino and the Southern Oscillation are different manifesta- 
tions of a single, large-scale, dynamic, air-sea, climatic feedback system (olten called "ENSO ). the 
inferred coupling has been somewhat less than perlect. For example, Deser and Wallace [52), based on 
the 1925 to 1986 interval, have found several instances when El Nino preceded the onsets of positive SLP 
at Darwin, Australia, and a few instances when El Nino followed the onsets ol positive SLP (including 
the largest El Nino of the century: the 1982-83 event). More peculiarly, they also found instances when 
El Nino and negative sw'ings in the SOI occurred separately. Thus, they concluded that El Nino and the 
Southern Oscillation are “more loosely coupled than previous studies would imply and that there are 
strong local controls on the climate ot the eastern equatorial Pacilic that sometimes transcend the 
influence of the Southern Oscillation." 

Models for predicting/hindcasting E.NSO are deterministic in nature; i.e., a change in an 
important oceanographic/atmospheric parameter (e.g., SST. SLP, SOL outgoing longwave radiation, 
surface wind stress, sea surface salinity, sea level height, the depth ol the thermocline, etc.) is either 
randomly or stochastically forced to bring about conditions lavoring the development ol an El Nino. 



While sometimes the forcing is introduced from the extratropies (via any one of a number of teleconnec- 
tion proeesscs), the loreing more often is introduced from entirely within the tropical Pacific region 
alone, where coupled interactions between the atmosphere and ocean are strongest [II, 16. 27-29, 35-38. 

66,70-71,80,87-88,98-99, 124-125, 153, 156-158, 174, 180, 183,207,216,248,260,275 288-289 
292, 317. 321-323], 

Because of the near simultaneity of two of the largest climatic events of the century, the early 
1982 eruptions of El Chichon in the northern tropical region of Mexico and the abrupt late- 1982 appear- 
ance of HI Nino (which followed El Chichon by about two seasons), much speculation has been raised 
that an association exists between them. In particular, because major volcanic eruptions are known to 
change the optical properties ol the Earth s atmosphere (especially the stratosphere), it has been sug- 
gested that alteration ol the atmospheric circulation occurs, leading either to conditions favorable for the 
initiation of El Nino (thus acting as the “trigger”) or to an intensification of an El Nino that is already in 
progress (thus modulating the severity of the El Nino) |94-96, 110. 196, 202. 229, 252, 257 1 

Some ol the statistical aspects of El Nino occurrences are examined here. In particular, the statis- 
tical properties ol the distribution of elapsed time between successive onsets of El Nino and the possible 
associational aspects ol El Nino occurrences, with respect to both volcanic and solar/geomagnetic 
activity, are examined. Reexamination of the association between volcanic activity and El Nino is 
prudent because recent reports have presented conflicting results [94-95, 196, 202 [ . Similarly, examina- 
tion of the association between solar/geomagnetic activity and El Nino, as yet to be explored, may yield 
results that are pertinent to El Nino prediction (in a “general” sense, not necessarily a “causal" sense). 
Finally, some remarks are given which may be of use for predicting the occurrence of the “next” 
anticipated Pd Nino. 


II. DATA 


Listings ol El Nino events can be found in Quinn et al. 1 221-222). Rasmusson [231 1, van Loon 
and Shea [278), Deser and Wallace |52|. and Ropelewski and Jones [ 247 [ , with the Quinn et al. [221 1 
listing being the most complete. The Quinn et al. [221 1 listing catalogs 79 El Nino events of near 
moderate-to-very strong intensity that appear to have occurred off “the west coast region of northern 
South America and its adjacent Pacific Ocean waters" during the interval 1525 to 1987, with events of 
moderate intensity (32) extending back to 1806. (No weak events are cataloged.) The Quinn et al. [221] 
listing represents a revision ol an earlier listing (Quinn et al. 1 222 ] ) and describes the very-strong, strong, 
and moderate events in the following manner: "The very strong events show extreme amounts of rainfall, 
Hood waters, and destruction, and coastal SST’s usually reach values 7° to 12 °C above normal during 
some months ol the southern hemisphere summer and fall seasons. The strong events, in addition to 
showing large amounts of rainfall and coastal Hooding and significant reports of destruction, exhibit 
coastal SST's in the 3° to 5 °C above normal range during several months of the southern hemisphere 
summer and fall seasons. The moderate events in addition to showing above normal rainfall, some flood- 
ing, and small amounts ol destruction, generally show coastal SST s in the 2° to 3 °C above normal range 
lor several months during the southern hemisphere summer and fall seasons. In all three categories the 
effects on coastal fisheries are highly damaging.” 
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Extensive listings of volcanic eruptions can he found in Lamb 1 I54|. Simkin et ^ 

and Sell 11941 Handler |94-961, and Bradley ( 1 9 1 . Major volcanic eruptions are usually desuibe *. 
h , nus Veil Index (DVD >?-»() or a Volcanic Explositivity Index (VEI) 5*4. As an example, 
^ -*15 (Tambora). ,he "lar g es, and deadbes, vnicanic e rap on 

; U- ,„ rv " (Slathers P67I) and had a DV1 of about 3.000. while the great eruption ot 1883 
(Krakatau) had a DVI of about lW according to Newhall and Self I '94| the ™ abated with these 

r-- 

eruptions, because of deposition on ice sheets in Greenland and Antarctica, ltkew.se. have been used to 
indicate the presence of major volcanic eruptions 151. 92-93, 160, 168]. as have estimates o atrnosp eric 
optical depth (Pollack et ah |215|). In the listings of Lamb 1 1541 and Newhall and Sell 1 194] the interva 
examined was from .500. Also, it should be noted that the volcanic record may be incomplete lor VEI 
=s4 for the period prior to about 1950 or 1960 1 194, 228 1. 

Sunspot (Zurich/ International sunspot number) and geomagnetic (the aa index) data ait ichably 
known backed about .848 and .868, respectively |61. .66-167, 176. 205 299 , while odaer solar/ 
geomagnetic parameters are of a more recent vintage. For example, the Ap index is known back U 
To. 7 -cm solar radio flux back to 1947. and satellite measurements ol the solar constant only back 
(Nimbus 7, N-7) or 1980 (Solar Maximum Mission, SMM). 

III. RESULTS 


Table 1 lists 39 El Nino events, from the aforementioned Quinn et al. U H listing, that oc curie 
between 1840 to 1987 and that were of nearly moderate-to-very strong intensity (strength. so i en i- 
fied in Table 1 are the particular years associated with each El Nino event an notes re erring ( 
from other studies (e.g , Rasmusson [231]; van Loon and Shea 12781; Ropelewsk. and Jones 1247 
These 39 events include 15 “strong” (S.S4, VS) and 24 “moderate” (W/M. M, M + > events^ A gencu 
note and a table legend appear at the bottom of Table 1. It is interesting to note that ol the _8 Warm 
Events cataloged by van Loon and Shea |278|, 24 are also found in Table 1 . The only ey ents appearing in 
van Loon and Shea that do not appear in Table 1 are for the years 1904. 1913, 1963 
the vears 1963 and 1969 were El Nino years, but the El N.ho events were ol weaker strength. Likewise, 
exc huhng event* 39 which does not appear in Rasmusson 1 23 1 1, of the remaining 38 events, r 1 appear in 
Rasmusson ’s listing, except events 5 (1860), 10 (1874). and 20 (.907); ,.e.. 35 out o 38 o the 
regionally determined El Nino events were also on a list ol El Nino events that were c ear y o g« 
proportions. Thus, while Quinn et al. [221 1 identified their El Nino events, trom which fable 
an extraction, primarily from local/regional effects, the events are truly global in extent. 

Table 2 lists 94 known volcanic eruptions that occurred between 1835 and 1986 that either had 
DVI 5? 50 VEI 2=4 or an observed dust cloud of column height 2*15 km and that may have had signifi- 
cant stratospheric effects. A column height of >15 km was selected as a criterion because in the tropics 
the tropopause (the dividing region between the upper troposphere and the lower stratosphere) averages 
above 15 km, while at mid and high latitudes, the height of the tropopause is below 15 km. The list ol 
volcanic eruptions was adapted primarily from lists tabulated by Lamb 1 154], Newhall and Seh |194] 
and Handler |94], with additions also from Eos (appearing in the column entitled "Ceoph^a Events : ) 
and from Toon and Pollack [272] (cf. Baldwin et al. ]9]; Pollack et al. ]215]) lor the period 1968 to 
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TABLE 1. 


SUMMARY OF EL NINO EVENTS: 1840-1987 
(From Quinn et al. [221]) 


Event No. (i) 


Year 


2 

1850 

M 

3 

1854 

W/M 

4 

1857-58 

M + 

5 

1860 

M 

6 

1864 

S 

7 

1866 

M 

8 

1867-68 

M 

9 

1871 

S + 

10 

1874 

M 

1 1 

1877-78 

VS 

12 

1880 

M 

13 

1884 

S + 

14 

1887-89 

W/M 

15 

1891 

VS 

16 

1 896-97 

M + 

17 

1 899- 1 900 

S 

18 

1902 

M + 

19 

1905 

W/M 

20 

1907 

M 

21 

191 1-12 

S 

22 

1914 

M + 

23 

1917 

S 

24 

1918-19 

W/M 

25 

1923 

M 

26 

1925-26 

VS 

27 

1930-31 

W/M 

28 

1932 

S 

29 

1939 

M ~f~ 


Strength Notes 

R: 1844-46; 2, 4, 3 
R: 1850; 2 
R: 1854; nc 
R: 1857; 2 
Not Listed 
R: 1864; 4; WE 
R: 1866; nc 
R: 1868; 3; WE (1868) 

R: 1871; 3 
Not Listed 

R: 1877-78: 4, 4; WE (1877) 

R: 1880; 3; WE 
R: 1884-85; 4,3; WE 

R: 1887-89; 3, 3, 1; WE (1888); Low (1888) 
R: 1891; 4; WE 

R: 1896; 3; WE (1896); Low (1896) 

R: 1899-1900; 4, 3; WE (1899); Low (1900) 
R: 1902; 3; WE; Low (1902) 

R: 1905; 3; Low (1905) 

Not Listed 

R: 191 1-12; 4, 3; WE (1911); Low (1911) 

R: 1914; 3 

R: 1917; 2 

R: 1918-19; 4, 3; WE (1918); Low (1918-19) 
R: 1923; 2; WE; Low (1923) 

R: 1925-26; 4, 4; WE (1925); Low (1925-26) 
R: 1929-30; 3, 3; WE (1930) 

R: 1932; 2; WE 
R. 1939; 3; WE; Low (1939) 
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TABLE I. (Concluded) 


Event No. (i) 

Year 

Strength 


Notes 

30 

1940-41 

S 

R: 

1940-41; 2, 4; Low (1939-41) 

31 

1943 

M + 

R: 

1943-44; 2, 2 

32 

1951 

W/M 

R: 

1951; 2; WE; Low (1951) 

33 

1953 

M + 

R: 

1953; 3; WE; Low (1953) 

34 

1957-58 

S 

R: 

1957-58; 4, 4; WE (1957); Low (1957) 

35 

1965 

M + 

R: 

1965; 3; WE; Low (1965) 

36 

1972-73 

S 

R: 

1972-73; 4, 4; WE (1972); Low (1972) 

37 

1976 

M 

R: 

1975-76; 1, 3; WE 

38 

1982-83 

VS 

R: 

1982-83; 4, 4; WE (1982); Low (1982-83) 

39 

1987 

M 

Not Listed 


Legend: W/M - Nearly moderate 

M - Moderate 
M + - Nearly strong 
S - Strong 

S + - Nearly very strong 
VS - Very strong 
R - Rasmusson [ 23 1 ] 

1 - Very weak 

2 - Weak 

3 - Moderate 

4 - Strong 

nc - Not classified 

WE - Warm event (episode) 

NOTE: The Rasmusson [231 1 listing, adapted from Quinn et al. [ 222 ] , also identifies El Nino events in 
1852 (2), 1855 (2), 1862 (nc), 1873 (2). 1875 (1), 1946 (1), 1948 (1), 1963 (1), and 1969(2). 
Of these, the events of 1963 and 1969 have been classified by van Loon and Shea [ 278 [ as warm 
events; the years 1904 and 1913 were also classified as warm events, although no counterpart 
El Nino event appears in either the Quinn et al. [221 1 or Rasmusson [231] lists for these two 
years. An event marked “Low" means that during the given year the Tahiti-Darwin Southern 
Oscillation Index remained in the lower 25 percent of the distribution for five months or longer 
(Ropelewski and Jones [247]); other “Low” years include 1913, 1946, 1963, 1969, and 1977. 
The term “Not Listed” means that the event does not appear in Rasmusson [231]. One should 
also note that on occasion there is a difference in the strength (intensity) classification in the 
various lists of El Nino. 
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TABLE 2. SUMMARY OF VOLCANIC ERUPTIONS WITH POSSIBLY SIGNIFICANT 

STRATOSPHERIC EFFECTS: 1835-1986 

(Adapted from Lamb [154], Newhall and Self [194], Handler [94], and descriptions reported in 
Eos (i.e. , dust clouds at 5* 1 5 km; Note: L-Lamb, NS - Newhall and Self, H - Handler, and E 
- Eos). Events noted as NS after 1869 also appear in Handler [94-95], and those events with 
possible well-defined signatures in either Greenlandic and/or Antarctic ice core deposits are 
denoted with an asterisk, taken from Hammer [92], Hammer et al. [93], and Legrand and 
Delmas [160].) 


Year 

Volcano 

Location (Degrees) 

Notes 

1835 

Coseguina 

13 N, 88 W 

NS; VEI = 5 (L; DVI = 4000); * 

1845 

Hekla 

64 N, 20 W 

NS; VEI = 4 (L; DVI = 800); * 

1846 

Am argu ra 

18 S, 174 W 

L; DVI = 1000; * 

1849 

Puraee 

2 N, 76 W 

NS; VEI = 4 

1852 

Etna 

38 N, 15 E 

L; DVI = 100 (50) 

1852 

Manna Loa 

20 N, 156 W 

L; DVI = 250 

1852 

Gunung Api 

4 S, 130 E 

L; DVI = 200 

1853 

Chikurachki 

50 N, 155 E 

NS; VEI = 4 

1854 

Sheveluch 

57 N, 162 E 

NS; VEI = 5 (L; DVI = ?) 

1855-56 

Cotopaxi 

1 S, 78 W 

E; DVI = 700; * 

1856 

Awu 

4 N, 126 E 

L; DVI = 70 

1856 

Sungihe 

4 N, 126 E 

E; DVI = 500? 

1856 

Komagatake 

42 N, 141 E 

NS; VEI = 4 (E; DVI <10) 

1861 

Makjan 

1 N, 128 E 

L; DVI = 800; * 

1869 

Puraee 

2 N, 76 W 

NS; VEI = 4 

1869 

Ceberuco 

21 N, 105 W 

E; DVI = 70 (50) 

1872 

Vesuvius 

41 N, 14 E 

E; DVI = 150 (70) 

1872 

Merapi 

8 S, 1 10 E 

L; DVI = 80 

1872 

Si nark a 

50 N, 155 E 

NS; VEI = 4 

1873 

Grimsvotn 

64 N, 17 W 

NS; VEI = 4 (L; DVI = 30) 

1875 

Askja 

65 N, 17 W 

NS; VEI = 5 (L; DVI = 1000); * 

1876 

Suwanose-Zima 

30 N, 17 W 

NS; VEI = 4 

1877 

Cotopaxi 

1 S, 78 W 

NS; VEI = 4 (L; DVI = 50) 

1881 

Nasu 

37 N, 140 E 

NS; VEI = 4 

1883 

Krakatau 

6 S, 105 E 

NS; VEI = 6 (L; DVI = 1000); * 

1883 

Augustine 

59 N, 153 W 

NS; VEI = 4 (L; DVI = 70) 

1885 

Falcon Island 

20 S, 175 W 

L; DVI = 300 

1886 

Tarawera 

38 S, 3 77 E 

NS; VEI = 5 (L; DVI = 800);* 

1886 

Tungurahua 

1 S, 78 W 

NS; VEI = 4 

1886 

Niafu 

16 S, 176 W 

L; DVI = 300; * 

1888 

Bandai san 

38 N, 140 E 

NS; VEI = 4 (L; DVI = 500);* 

1888 

Ritter Island 

6 S, 148 E 

L; DVI = 250 


TABLE 2. (Continued) 


Volcano 

Suwanose-Zima 

Bogoslof 

Awu 

Thompson Island 
Una Una 
Dona Juana 
Mont Pelee 

Soufriere (St. Vincent) 
Santa Maria 
Thordarhyma 
Bogoslol 

Shtyubelya 
Sopka (Ksudach) 


Location (Degrees) 

30 N, 130 E 
54 N, 168 W 
4 N, 126 E 
54 S, 5 E 
0 S. 121 E 
2 N, 77 W 
15 N. 61 W 

13 N, 61 W 
15 N, 92 W 
64 N, 18 W 
54 N, 168 W 
52 N, 158 E 


1909 

Tarumai 

43 N, 141 E 

191 1 

Taal 

14 N, 121 E 

1912 

Katmai (Novarupta) 

58 N, 155 W 

1913 

Colima 

19 N. 104 W 

1914 

Sakurashima 

(Sakurazima) 

32 N, 131 E 

1917 

Agrigan 

19 N, 146 E 

1918 

Tungurahua 

1 S, 78 W 

1918 

Katla 

64 N, 19 W 

1919 

Man am 

4 S, 145 E 

1921 

Andes (Puyehuej 

41 S. 72 W 

1924 

Raikoke 

48 N, 153 E 

1929 

Komaga-Take 

42 N, 141 E 

1931 

Kluchev (Kliuchevskoi) 

56 N. 161 E 

1932 

Fuego 

14 N, 91 W 

1932 

Quizopu (Cerro Azul) 

36 S. 71 W 

1937 

Rabual 

5 S, 152 E 

1945 

Kliuchevskoi 

56 N. 161 E 

1946 

Sarychev 

48 N, 153 E 

1947 

Hekla 

64 N. 20 W 

1951 

Mt. Lamington 

9 S, 148 E 

1951 

Ambryn 

16 S, 158 E 

1952 

Bagana 

6 S, 155 E 

1953 

Mt. Spun* 

61 N, 153 W 

1955 

Nilahue 

40 S, 72 E 

1956 

Bezymjannaja 

(Bezymianny) 

56 N, 161 E 


NS; VEI = 4 
L; DVI = 1 70 (50) 

L; DVI = 100 

L; DVI = 1300 (400); *; Island destroyed 
L; DVI = 140 
NS; VEI = 4 

NS; VEI = 4 (L; DVI = 100); 

Two eruptions of VEI = 4 
NS; VEI = 4 (L; DVI = 300); * 

NS; VEI = 6 (L; DVI = 600); * 

NS; VEI = 4 
L; DVI = 70 (20) 

NS; VEI = 5 (L; DVI = 500); * 


NS; VEI = 4 

NS; VEI = 4 (L; DVI = 30) 

NS; VEI = 6 (L; DVI = 500); * 

NS; VEI = 4? (L; DVI = 10) 

NS; VEI = 4 (L; DVI = 40) 

NS; VEI = 4 
NS; VEI - 4 
NS; VEI - 4 
NS; VEI == 4 

NS; VEI == 4 (L; DVI = 200) 

NS; VEI ■ = 4 
NS; VEI = 4 

NS; VEI = 4 (L; DVI = 5) 

NS; VEI = 4 

NS; VEI = 5 (H; 6) (L; DVI = 70) 
NS; VEI = 4? 

NS; VEI = 4 
NS; VEI = 4 

NS; VEI = 4 (L; DVI = 70); * 

Dust cloud to 27 km 

NS; VEI = 4 (L; DVI = 20) 

NS; VEI - 4 
NS; VEI = 4 

NS; VEI = 4 (L; DVI = 7); 

Dust cloud to 23 km 

NS; VEI = 4 

NS; VEI = 5 (L; DVI = 30) 

Dust cloud to 45 km 



TABLE 2. (Concluded) 



1960 

Punliagudo 

42 S, 72 W 

1963 

Gunung Agung (Agung) 

8 S, 1 16 E 

1964 

Sheveluch 

57 N, 162 E 

1965 

Taal 

14 N, 121 E 

1966 

Kelut 

8 S, 1 12 E 

1966 

Awu 

4 N, 126 E 

1966 

Oldoinyo Lengai 

3 S, 36 E 

1968 

Femandina 

0 S, 92 W 

1973 

Tialia 

44 N, 146 E 


1974 

1975 

1976 

1977 
1979 

1979 

1980 


Fuego 

Plosky Tolbachik 

Augustine 

Bezymianny 

Bezymianny 

Soufriere 

Mt. St. Helens 


14 N, 9! W 
56 N, 160 E 
59 N, 153 W 
56 N, 161 E 
56 N, 161 E 
13 N, 61 W 
46 N, 122 W 


1980 

Hekla 

64 N, 20 W 

1980 

Gareloi 

52 N, 179 W 

1981 

Alaid 

51 N, 156 E 

1982 

Mystery 


1982 

El Chiehon 

17 N, 93 W 

1983 

Mystery 


1984 

Home Reef 

19 S, 175 W 

1985 

Nevado del Ruiz 

5 N, 75 W 

1986 

Pavlof 

55 N, 162 W 


L; DVI = 100 (50) 

Several other vents (unnamed); 

Given latitude/longutude is average 

NS; VE[ = 4 (L; DVI = 800); * 

NS; VEI = 4 

NS; VEI = 4 (L; DVI = 10-15) 

NS; VEI = 4 

NS; VEI = 4 (L; DVI = 50-100) 

Dust cloud to 18 km 

NS; VEI = 4 

NS; VEI = 4 (L; DVI = 50-100) 

NS; VEI = 4 
NS; VEI = 4 
NS; VEI = 4 
NS; VEI = 4 
E; Dust cloud to 15 km 
NS; VEI = 4 
E; Dust cloud to 18 km 
NS; VEI = 5; 

E; Dust cloud to 3=19 km 
E; Dust cloud to 15 km 
E; Dust cloud to 10.5 km; S0 4 at 19.2 km 
E; Dust cloud to 15 km 

E; Dust at 10-20 km (January-March); 

Several eruptions possibly associated 
with cloud 

H; VEI = 4; E; Dust cloud to 16.8 km 
(possibly to 26 km) 

E; Dust at 18-19 km (February-April) 

E; Dust cloud to 15 km 

E; Dust cloud to 10.5 km (possibly to 25 km) 
E; Dust cloud to 16 km 


spheric" TOl^lT l ^IJ^^S l ^an^ , ^l ha ^Sr«srd-ff Snir f ;ant T™* ° f S ‘ rat °' 
(Handler ,94.95,, E tha , 'J 

stratosphenc aerosols (1973 - Tiatia; .975 - Plosky 

ezymianny) while 10 were added (August 1928 - Paluweh, 8 S 121 W 1929 

Reventador 0 S, 77 W; September 1971 - Fuego, .4 N 91 W- April 1979 

Soufriere, 13 N 61 W- November iq 7 q c- kt , ’ 1 w » A P nI iy '9 - 

it. ’ ’ vv, iNovember 1 979 - Sierra Negra, 1 S 91 W Oomhor jocn 

oeen me source ot the Mystery event of early 1982. J 



The year of the eruption, its name and location, and notes regarding each of the entries are given. An 
asterisk (*) denotes those events which also have a reported signature in the ice core deposits of either 
Greenland or Antarctica. Of the 94 events, 92 can be associated with specific volcanic eruptions of 
known location; two are denoted "Mystery." (The mystery cloud of 1982 may be associated with 
Nyamuragia, located at 1 S, 29 E, that erupted in December 1981; Handler [96].) Excluding the two 
mystery aerosols, 45 occurred within the tropics (within 23.5 degrees of the equator, a region containing 
about 40 percent of the Earth's surface area, predominantly water, as well as the bulk of the Earth's 
atmosphere). In particular, 23 occurred in the northern tropics, 22 occurred in the southern tropics, 41 
occurred north of 23.5 N (the bulk of which occurred north of 45 N), and 6 occurred south of 23.5 S (the 
bulk of which occurred between 23.5 S and 45 S). Also, of the 94 eruptions, 17 have been clearly identi- 
fied in ice core deposits (perhaps, 18; evidence has been presented by Hofmann et al. [ 1 14] that the El 
Chichon stratospheric aerosol was detected in Antarctica). 

Plotted in Figure 1 are the annual averages of sunspot number, l().7-cm solar radio flux, the aa 
and Ap geomagnetic indices, and solar irradiances (SMM and N-7) for data available periods during the 
interval 1840 to 1987. The long-term average for each of these parameters is drawn as the horizontal line. 
As an exammple, 54.3 represents the average annual sunspot number for the interval 1 848 to 1987, 1 29.5 
represents the average annual 10.7-cm solar radio flux value for the interval 1947 to 1987, and so on. 
Individual sunspot cycles are identified al the bottom of Figure I . Across the top of Figure 1 are the 39 El 
Nino events given in Table 1 and the 94 possibly significant volcanic eruptions given in Table 2, where 
the shaded triangles denote “strong" El Nino events and tropical volcanic eruptions, respectively. The 
term “VS” means that the El Nino was very strong. (The mystery aerosol of 1982 has been shaded based 
on the remarks previously given; Handler [ 96 1 . ) 

In the next section, discussion will center on Figure 1 (and the supporting Tables I and 2). It is 
upon this set of comparative results that the bulk of the remaining text is concerned. 


IV. DISCUSSION 


A. Statistical Aspects of El Nino Onsets 

From Figure 1 and Table I , one notes that 38 El Nino events of nearly moderate to very strong 
intensity occurred during the 140-year history of modern era sunspot observations (1848 to present). Of 
these 38 events, 24 have been of moderate (W/M, M, M + ) intensity and 14 of strong (S, S + , VS) 
intensity. Thus, on average, El Nino occurs at the rate of one event per 3.7 years or, incorporating 
intensity, at the rates of one moderate event per 5.8 years and one strong event per 10 years. It may also 
be noteworthy that every decade between the I840's and the I980’s, except the 1960's, had a moderate 
to stronger El Nino that spanned two consecutive years (e.g., 191 1-12, 1925-26, etc.; see Table 1) and 
that there has never been a single instance of back-to-back strong El Nino events. 

In Figure 2 (top-left), the elapsed time (in years) between successive onsets of El Nino is plotted 
as a function of event number (see Table 1 ). As noted above, the mean interval between onsets of El Nino 
is about 3.7 years, identified in Figure 2 as the horizontal line running parallel to the event number axis. 
The sample of 38 elapsed times has a standard deviation of about 1 .8 years and a mode and median of 3 
years. A runs test (Lapin [155], p. 624) suggests that the sample of elapsed times is randomly distributed. 
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Figure 1. Annual counts of El Nino in relation to annual counts of possibly significant volcanic 
eruptions and to annual averages of sunspot number, 10.7-cm solar radio flux, aa and Ap 
geomagnetic indices, and solar irradiance (Nimbus 7 and Solar Maximum Mission) for 
the interval 1840-1987. Very strong (VS) El Nino events are denoted, as are several of 
the major volcanic eruptions (by name). Long-term averages are identified. See text for 
other details. 
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Figure 


Elapsed time (in years) between successive dls^r'ibutloTof 

and the distribution of El Nino strengths (bottom-right). 
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STRENGTH 


although greater variation is clearly evident from about event number 23 (191 7) than before it Subdivid 

i £ s,ze f samp,e v 2 to 20 and 21 *»• 

g t up (_ to .0) has a mean elapsed time of 3.3 years and a standard deviation of I 1 years while the 

T mng 8r ° UP l 2 ' l ° 39) has 3 mean da P sed time of ' 4 - 2 years and a standard deviation of 
■ ^f° esis tes bng, based on the t statistic for independent samples (Lapin [1551 p 486) and 

assummg each sample to be normally distributed, suggests that the two samples are not signiHcantly dif 
ent, implying that the two groups can be combined into a single grouping and that a dnoi^ ^ ■ a 

rr deVi r i0n Can be USed t0 deSCribe the ° Vera " set of ela P sed times. (It Should be noted. however, 
neak- at e t0n ^ dmmg events 2 t0 20 has elapsed times spanning from I to 6 years, including a stron" 

mrlM f y ! arS ’ W a 6 C ° ntainin§ events 21 to has elapsed times spanning 1 to § 8 years^ 

ncludmg a less striking peak at 3 years. Thus, while the former sample has a distribution that looks like 

hat ofTh f d° n, K C attCr Sample ’ S distribution does instead, it has a distribution more like 
th m r" 1 ° rm dlStr i ' bUt,0n ‘ A|S °’ beCause greater fiance exists in the latter sample as compared to 

than'the ° f even, » 21 is larger 

variable climate may have accused during ,he present centu^, IZ7 

o e wonders whether the change is tnerely a natural ept.sodic swing in dimate or is i, an a^hrop g „ e 


Also plotted in Figure 2 (top-right) is the overall distribution of elapsed times making no distinr 
"°n betw “" *«»«? ■» 20 and 21 39. Because the mode and median are smaTr Z 2: meat o„e‘ 

perceives that the distribution is positively (or rightward) skewed. Application o( the Kolmogorov 

3™ w and'^ot eVla, '° n 'IT E °° dneSS ° r fU <Lapi " f'551. p. 639), using an assumed value of 

respectively, lor the mean and standard deviation of the population of elapsed 
tmes. suggests that the distribution of elapsed times is normally distributed; on the other hand a p Inula 

„«tt 0r4 „ ,, TT'u 2 0 yearS ' reS|KC,iVdy ' ^ *"a, Jhe distrZ on 

, P ' llieS md l not be normally distributed. (The distribution of elapsed times for the entire 

:r„S’ If’ 7 F,gUre looks "0™ally distributed; based on a sample distribution of 38 

elapsed times, having a sample mean of 3.7 years and a sample standard deviation of I 8 velr, one 
computes the population mean and standard deviation to be |x = 3.7 + 0 6 and o = | 8 +06/0 ^ 
respectively, at the 95 percent level of confidence.) +0.6/-0.3, 

„ n Jl !l abOVe t analy f is L SUggests that El Nino events of moderate to stronger strength do not occur on 
" Ua as ‘ S ' nstead ’ the Y usually occur about once every 3 years (one third of the entire sample- 
two-thirds of the entire sample have elapsed times of 2 to 4 years) and only occasionally do they occur in 

pi °, aC y ^ JrS e a P seti h 1116 equals 1 year which is observed in only about 8 percent of the sample) 

f M ° ! ha ' CXCeed 5 yearS have ^ —pi a« h o e f ir an 

as eviaent early on, with 5 of the 6 events occurring in the last 50 years Because the onsets of FI Mi™ 

usua ly occur alter the Earth s closes, approach the Sun when inflation is grelleM ofsup^s 

msolal t, on should be a major lactor in the genesis of El Nino; however, becaufe of the lack of an annual 

signal ,n the elapsed times between successive El Nino of moderate to stronger intensity one surmises 

dur nHo th n 3 h ne ,Ully CXpla,n thC observed variation. (El Nino events usually have their onsets 

,• .V Uther " hemispheric summer/fall, corresponding to December-May; a comparison of the events 

tsted in Table I to those identified in Deserand Wallace |52| forthe interval 1978 to 1986 ,ho ,u 

tan 80 percent of the El NiDo events had their onsets doling ”Z 

greatest number during southern summer and with no events during southern winter-June-August. 
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It should also be noted that an annual signal occurs in the 10.7-cm solar radio flux, which is a conse- 
quence of the Earth's orbital motion, being higher, on average, during January than during July, and a 
semiannual signal occurs in geomagnetic activity, which may be a consequence of the tilt of the Sun’s 
rotational axis relative to the equator, with enhanced geomagnetic activity in March and September, 
when more of the Sun's polar regions face the Earth, than in June and December.) 

Also plotted in Figure 2 (bottom-left and bottom-right, respectively) is the time variation of the 
relative strengths of the moderate to stronger El Nino events and their distribution. Although the 
strengths appear to vary randomly, moderate El Nino events occur more often than strong El Nino events 
and the time series of El Nino strengths has, as yet, never included a strong-strong occurrence. Thus, if 
one were forecasting the relative strength for the “next” El Nino, one probably would forecast a moderate 
event, especially if the last known event was a strong El Nino; on the other hand, if the last event was a 
moderate event, one could only say that the next anticipated event will be of moderate to stronger 
strength (ignoring the weaker El Nino events). 


B. El Nino Onsets Versus Volcanic Activity 

More than 200 years ago, Benjamin Franklin [ 72 ) noted a coolness during the summer months of 
1783 that was associated with a “constant fog” over Europe and North America. According to Franklin, 
surtace temperature decreased because some of the incident sunlight could not reach the Earth's surface, 
owing to the presence of the fog. Of several possible causes for the fog. Franklin speculated that volcanic 
activity in Iceland might be responsible for the meteorological effects |23, 50, 215, 254, 261], 

Later, Humphreys 1 1 2 1 J showed that preferential scattering of sunlight by volcanic aerosols could 
alter the Earth’s albedo, thereby effecting climatic changes (cf. Humphreys [122]), and Kimball 1 139- 
1 40 j showed that reduction in the amount of direct sunlight at the Earth’s surface by as much as 20 to 30 
percent of that ordinarily observed was seen after certain volcanic eruptions [190, 290-291 ]. In 1961 a 
persistent, world-wide, stratospheric aerosol layer, composed mainly of sulfate particles (and sulfuric 
acid droplets), was detected by high-altitude balloon and aircraft sampling experiments |22, 129-130]. 
Results ol monitoring this layer have since shown that direct volcanic injection occasionally takes place, 
having both immediate short-term and severe longer term effects [79, 181], As an example of a possible 
climatic ettect related to volcanic eruptions, Kondo [ 143) has noted that “most of the poor rice harvests 
caused by unusual cool summers and leading famine conditions in the Tohoku district (the northeastern 
part ol Japan) took place just after the great volcanic eruptions” during the last 300 years (cf. 1 19, 94-96, 
163, 188, 196, 202 ] ) . It is perceived, then, that volcanic eruptions, especially violent ones in the tropics, 
because of their globally perturbing aerosols (Lamb [ 1 54 1 ) , alter the Earth’s atmospheric circulation, 
thereby, influencing both local and world-wide climatology [5-6, 30-32, 46-47, 53, 58-59, 76, 82-83, 
97, 100, 111, 113, 118, 120, 123, 127, 133-134. 150, 152, 159, 164, 170-172, 175, 182, 184-185. 187, 
191, 193, 197, 203-204, 214, 223, 229-230, 235, 239-243. 245, 255-256, 258, 266, 269, 281-283, 
312], 


Recently, Handler [94-95] examined the relationship between El Nino occurrences and volcanic 
eruptions, to determine whether or not the 1982-83 El Nino was related to the March-April 1982 
eruptions of El Chichon. In his analysis, he composited (using the superposed epoch analysis technique; 
cf. Haurwitz and Brier f 1 03 1 ) separately and seasonally the SST’s for 1 1 low-latitude (<20 degrees) and 
20 high-latitude (>20 degrees) eruptions occurring between 1868 and 1980, where the SST’s were for 
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the eastern equatorial Pacific bounded by 0 to 10°S and 90 to I80°W. His results suggest that, when 
major stratospheric aerosols occur at low latitudes, the SST’s ot the eastern equatorial Pacific tend to be 
warmer than normal for the three seasons following the eruptions and that, when major stratospheric 
aerosols occur at high latitudes, the SST’s tend to be cooler than normal for about four or five seasons 
following the eruptions, concluding that the association between SST anomalies and volcanic eruptions is 
statistically significant at the 95 percent level ot conlidence. 

While Parker [202], using a different data set, has confirmed Handler's findings, Nicholls [ 196] 
notes that, in composites (10 events) of Darwin SLP, another measure often employed to indicate the 
occurrence of El Nino, the date of eruption tends to be preceded by lower than average pressure and 
followed by higher than average pressure. He further found that “a strong linear upward trend in the 
composite pressure anomaly starts well before the date of eruptions and continues for several months" 
afterward, suggesting that ENSO events do not result from low-latitude volcanic eruptions (owing to 
violation of causality). Similar conflicting results can be found in Sear and Kelly [ 257 ] both tor 
individual and composited events based on northern hemisphere air temperature (cooling is seen to have 
begun some 3 months prior to the El Chichon eruption and 5 months prior to the dates of eruption based 
on a composite of five eruptions). 

To assess the associational aspects of El Nino and volcanic eruptions, from Figure 1 , one can 
compare the onsets of El Nino with the occurrences of volcanic eruptions. From the preceding, it a strong 
association exists between volcanic eruption in the tropics (cause) and the onset of El Nino (effect), then 
one expects volcanic eruptions to preferentially occur before (or simultaneously with) the appearances of 
El Nino, because the response should be essentially immediate with only slight delay (within three to five 
seasons; certainly, within 2 years). (Observations ot the decay ot the aerosol load following selected 
eruptions indicate that it has an e-folding time of about 1 year, implying that the original load has been 
reduced by about a factor ot 2.7, 1 year alter the eruption, by about 7 .4, 2 years after the eruption, and by 
about 20.1, 3 years following the eruption. Also volcanic materials have been observed to circum- 
navigate the globe in less than 1 month and to span from pole to pole in about I year, with the maximum 
opacity usually occurring within the hemisphere of origin and in less than 1 year following the eruption.) 


Figure 3 plots the results of a comparison of occurrences of volcanic eruption relative to El Nino 
onsets for tropical volcanic eruptions alone (bottom) and for all major volcanic eruptions (ignoring the 
latitudinal position of the eruptions; top). If El Nino events result purely from the eruptions ot major low- 
latitude (tropical) volcanoes that precede the onsets of El Nino by =£ 1 year (within seven seasons), one 
finds that only 16 out of 38 El Nino events (between 1848 and 1987) can possibly be "explained.” If one 
opens the window to include all major tropical eruptions within 2 years (within 1 1 seasons), then one 
finds that the supposed association can only account lor 22 out of 38 El Nino onsets (about 60 percent). 
Thus, nearly half (2*40 percent) of the El Nino onsets are not preceded (within 2 years) by a major 
volcanic eruption in the tropics. 

From Figure 3 (top), if one ignores the location of the volcanic eruption, one finds a stronger 
apparent association between them. For example, within 1 year, 26 out of 38 El Nino onsets were 
preceded by a major eruption and, within 2 years, 30 out of 38 El Nino onsets were preceded by a major 
eruption; only 8 out of 38 (about 20 percent) were not preceded by a major eruption. Thus, by ignoring 
latitude, one can attribute more of the El Nino events to the occurrences ot major volcanic eruptions; 
however, the association does not appear to be one-to-one. (It should be emphasized that this analysis 
was performed relative to El Nino onset.) 
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Figure 3. 



ELAPSED TIME (YR) 
BETWEEN VOLCANIC 
ERUPTION & EL NINO ONSET 


«,<„„> and latitude-independent (denoted “All;" top. major volcame erupt, ons. The 
ished lines refer to events occurring since 1950. 
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Limiting oneself to just those events that have occurred since 1950 (8 events, the dashed lines in 
igure 3), one finds that 6 out ol 8 El Nino onsets were preceded (within 2 years) by major volcanic 
eruptions in the tropics, or that the association is roughly in the same proportion as that of the larger 
r mpanson. gnoring latitude, one tinds 7 out of 8 El Nino onsets to be preceded (within 2 years actuallv 
within 1 year) by major volcanic eruptions. Thus, although the more reliable data set indicates a possible 
connect, on between E. Nino onsets and preceding major volcanic eruptions, the i„W S 
remains not one-to-one. Event 36 (Table 1 ), the strong E. Nino of .972 73, cannot be attribu ted to any 
known major volcanic eruption (s£3 years), regardless of latitude. Y 

htitn t T? b ' e - 3 n' StS . thep0Ssib,e ass °ciations that may exist between El Nino onsets and major low- 
tude (tropical) volcanic eruptions for elapsed times spanning 0 to 3 years. It identifies the particular El 
o event, its onset year and strength, and the name, location, VEI (if known; otherwise^ DV1 is 

ofFiiure^ ToVnIuT^ ^ 7r SSib !? associated volca "°- As previously noted in the discussion 

eruptions OHhese 1 0 T" T 8 7 Same years of the 0 «urrences of major, tropical volcanic 

eruptions. Of these 10, 8 occur during the first part of the year when El Nino is usually observed- two 

7 m thC makUlg thCSe associations less Kkely. These two events are numbers 17 
and 35. Instead of same year occurrence for both the El Nino and the eruption for these two events event 

7 may be better associated with the eruption of Una Una (0°S) in 1898, implying an elapsed time of 

. ye . dr ’ and ev f nt 35 ma y be better associated with the eruption of Gunung Agung (8°S) in March 1963 
implying an elapsed time of 2 years. ' 

Table 4 identifies all of the extratropical, southern hemispheric eruptions (6) that occurred nrior tn 
the onsets of El Nino. While it is less likely that these eruptions produced Ltosp^^^^ 

8 f ° a 7 eX 7 t(t C d d f e ls that extratro P ] cal eruptions tend to be most effective within their hemisphere 
DverTssiWh C 7 rUptl °" S Can be effective throughout the tropics, even globally; Lamb [154] 

lufthl ^ rcL Ti°ui haVC P roduced aerosoIs that ^re effective in the southern hemisphere. Combin’ 

L a es ~ anc ^ 4, one that, it El Nino can be triggered by preceding (within 1 vearO 
major, tropical and/or extratropical (southern hemisphere) volcanic eruptions then such an association 
can only account for about two-thirds of the observed El Nino (25 out of 38 events); one-third must be 
nggered by something else. As has been previously noted, the association appears stronger if one 

assodM ^ eW ' Se ’ by ,ncludin g smaller eruptions (VEI =S3), one also can infer a stronger 


C. El Nino Onsets Versus Solar/Geomagnetic Activity 

Two examples of solar-terrestrial relationships are the distribution of average air temperature 
from the equator to the pole (wanner near the equator as compared to cooler near the pdes) and the exist 
ence of aurorae the implication being that a changing Sun is, perhaps, a major piece of the puzzle 
regardmg weather/climate. While the Sun is recognized as the ultimate source of energy for the 
atmospheric-oceanic climatic system, because a full understanding of the step-by-step processed at work 

w.thin the system has not yet been realized, researchers remain cautious, even skeptS towards any 
claims as related to Sun-weather/climate. 1 rowaras any 
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Event Onset 


Strength 


Elapsed time (year) = 0 


17a 

18 


21 

23 

24 
28 
32b 

35c 

38 


1877 

1899 

1902 

191 1 

1917 

1918 
1932 
1951 

1965 

1982 


Volcano 


VS Cotopaxi 

S Doha Juana 

VI + Santa Maria 

Soufriere 
Mont Pelee 

S Taal 

S Agrigan 

W/M Tungurahua 

S Fuego 

W/M Ambryn 

Mi. Lamington 

M + Taal 

VS El Chichon 


Latitude 


1 S 

2 N 

15 N 

13 N 

15 N 

14 N 
19 N 
1 S 
14 N 

16 S 
9 S 

14 N 

17 N 


VE1 


4 

4 

6 

4 

4 

4 

4 

4 

4 

4 

4 

4 


a. Possibly better associated with 

b. Possibly better associated with 

c. Possibly better associated with 


Una Una (Lat. 0 S) which erupted in 1898. 
Mt. Lamington than Ambryn. 

Gunung Agung (Lat. 8 S) which erupted in 


Month 


June 

November 

May 

May 

May 

June 

April 

April 

January 

September 

January 

September 

April 


March 1963. 


Elapsed time (year) — 1 


2 

1850 

M 

4 

1857 

M + 

13 

1884 

S + 

14a 

1887 

W/M 

22 

1914 

M + 

33 

1953 

M + 

a. 

Tarawera (38 

S) erupted 


Purace 

2 N 

Sangihe 

4 N 

Awu 

4 N 

Cotopaxi 

1 S 

Krakatau 

6 S 

Niat'u 

16 S 

Tungurahua 

1 S 

Colima 

19 N 

Bagana 

6 S 

; 1886, having VE1 

= 5. 


4 December 

(DVI = 500?)March 
(DVI = 70) ? 

(DVI = 700)7 
6 August 


(DVI = 300) ? 

4 June 

4 June 

4 February 
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TABLE 3. (Concluded) 


Event Onset 


Strength 


Elapsed time (year) = 2 


Volcano Latitude 


3 

1854 

W/M 

9 

1871 

S + 

10 

1874 

M 

29 

1939 

M + 

37 

1976 

M 

39 

1987 

M 

a. 

Dust cloud possibly to 25 

Elapsed time (year) 

= 3 

6 

1864 

S 

12 

1880 

M 

15 

1891 

VS 

19 

1905 

W/M 

30 

1940 

S 


Gunung Api 

4 S 

Purace 

2 N 

Merapi 

8 S 

Rabual 

5 S 

Fuego 

14 N 

Ruiz 

5 N 


Makjan 

0 N 

Cotopaxi 

1 S 

Ritter Island 

6 S 

Santa Maria 

15 N 

Soufriere 

13 N 

Mont Pelee 

15 N 

Rabual 

5 S 


VEI Month 


(DVI = 200) ? 


4 October 

(DVI = 80) April 
4 May 

4 October 

a November 


(DVI — 800) December 
4 June 

(DVI = 250) ? 

6 May 

4 May 

4 May 

4 May 


TABLE 4. MID-TO-HIGH LATITUDE SOUTHERN HEMISPHERE VOLCANIC FRUPTiniue 

that possibly may be associated with el NISO events S 


Year Volcano 


Latitude 


VEI 


1886 Tarawera 38 S 

1895/96 Thompson Island 54 S 

1921 Puyehue 41 s 

1932 Quizopu 36 S 

1955 Nilahue 40 S 

I960 Puntiagudo 42 S 


(DVI = 1300) 
4 
6 
4 

(DVI = 100) 


Month 

Event 

Year 

Strength 

June 

14 

1887 

W/M 

9 

16 

1896 

M + 

December 

25 

1923 

M 

April 

28 

1932 

S 

July 

34 

1957 

S 

May 

35 

1965 

M + 
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One of the first in this century to suggest a connection between the Sun and terrestrial weather/ 
climate was Humphreys [120], who also related the production of terrestrial ozone to solar ultraviolet 
(UV) radiation and to auroral discharge. Variations over the solar cycle of the UV flux, the solar con- 
stant, and the geomagnetic activity (related to the occurrences of solar Hares, coronal holes, erupting 
prominences, and even the time of year) have been documented which have been associated, in par- 
ticular, with changes in the upper and middle atmospheres and, possibly, with changes even in the lower 
atmosphere, although it is this latter aspect that is often controversial 1 2-4, 10, 17, 21 , 24-26, 33, 42-45, 
49, 54-57, 73, 75, 77, 81, 84-85, 90, 102, 104-105, 107-109, 115, 126, 131-132, 141-142, 144, 149, 
151. 169, 177, 179, 186, 198-200, 205, 21 1-213, 217-219, 236-237, 250-251, 253, 259, 263-265, 268, 
270-271, 277, 279-280, 293-297, 312, 324], 

The term “solar activity cycle” (or solar cycle) is the general term employed to describe the 
multiplicity of solar-related phenomena that appear to vary cyclicly in a periodic or quasi-periodic 
manner, having a period of approximately I I (or 22) years. The oldest measure ol the solar cycle is the 
relative sunspot number first introduced by Rudolf Wolf in 1848 (cf. Waldmeier [285]; Kiepenheuer 
[137]; Wilson [ 299 1 ; McKinnon 1 176]). Other measures include, for example, the 10.7-em solar radio 
flux, sunspot area, photospheric line strengths, number of sunspot groups, number of sunspots, number 
of solar Hares, etc. In a larger sense, the solar cycle includes not only the sunspot cycle, but also 
variations of the solar irradiance (which is believed to vary over the sunspot cycle; e.g., Hudson 1 1 1 9 1 ; 
Livingston, Wallace, and White 1 161 1: Willson and Hudson [298]) and the geomagnetic cycle (since the 
magnetic disturbances that give rise to the fluctuations, especially as seen in the aa and Ap indices, are 
believed to be streams of charged particles of solar origin; cf. Garland 1 78 ] ; Feynman and Gu [ 68 ] ; 
Mavromichalaki, Vassilaki, and Marmatsouri |165|). Provided that statistically significant correlations 
exist between the various parameters and sunspot number, one can use the sunspot cycle as a proxy tor 
the solar cycle, thereby allowing a convenient way to assess statistical aspects of El Nino, especially 
those that might be suggested as being indicative of preferential behavior. 

Scatter plots of solar irradiance (from Nimbus 7 and Solar Maximum Mission), geomagnetic 
indices (aa and Ap), and the l().7-cm solar radio llux against sunspot number (using annual averages) are 
shown in Figure 4. In each case the depicted linear fit is found to be statistically significant, although 
considerable variation is seen in the geomagnetic indices. On the basis of these linear fits, one infers that, 
indeed, the various parameters are directly correlated with sunspot number. Thus, sunspot number (the 
sunspot cycle) can be used as a proxy for the solar/geomagnetic cycles and for the study of statistical 
associations between onsets of Fd Nino and the solar cycle. Statistically speaking, one infers that >80 
percent of the variation observed in the l().7-cm solar radio flux can be “explained” by the variation 
observed in sunspot number, while >90 percent of the observed variation in the solar irradiance as 
measured by the Solar Maximum Mission can be explained by the variation observed in sunspot number. 
The percentage of explained variation is about 20, about 30, and about 50 percent for the Ap index, the 
aa index, and the solar ii radiance as measured by the Nimbus 7, respectively. 

The results of superposed epoch analyses, using sunspot minimum occurrence as the epoch ol 
comparison, are illustrated in Figure 5 for Ap, aa, 10.7-cm solar radio llux, and sunspot number. As 
noted above, the geomagnetic indices, the 10.7-cm solar radio llux, and sunspot number all display a 
somewhat cyclic appearance, having a period of about 10 to 1 1 years and a maximum about midway 
through the cycle. More specifically, while sunspot number and 10.7-cm solar radio (lux have minima 
that coincide exactly, Ap and aa minima usually follow sunspot minimum by about 1 year. For Cycles 18 
and 19, maxima in sunspot number and 10.7-cm solar radio llux coincide exactly, while lor Cycles 20 
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Figure 4. Linear fits of 10.7-cm solar radio flux, aa, Ap, and solar irradiance (Nimbus 7 and Solar 
Maximum Mission) against sunspot number using annual averages of each. The linear 
regressions are shown as are the correlation coefficients (r). 
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Figure 5 
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Sane, nosed epoch analyses of sunspol number. 10.7-cm solar radio llux aa. and Ap 
b sS P ": L epoch of y sunspo, number min, mum (using annual averts . — 
cycle minima and maxima are identified lor each parameter. The as ensh < ^ 

rvrle n indicates that for both aa and Ap. a minimum occurred at the relative position 
shown' however, the lowest minimum occurred in 1980. near sunspot max, mum. uni, he 

any previous cycle. 
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and 2 1 the radio (lux maximum followed the sunspot number maximum by about 2 years Maxima in the 
o geomagnetic indices have always coincided with each other (as have their minima) and thev hwp 

I*n C°vcle f ° IIOWed their res P ective suns P ot maxima by two or more years with only two exceptions' 

the i I h- maX,ma S 7 Sp0t " Umber and the geomagnetic indices coincided and in Cycle 13 

he geomagnenc maximum preceded the sunspot maximum by 1 year. Sunspot maximum has always 

sol ar r!di Tn ^ by 4 ± 1 year for the most reliably known sunspot cycles, and the 10 7-cm 

solar radio flux maximum has always followed its minimum by 3 to 6 years. The difference in epochs of 

maximum for sunspot number and the 10. 7-cm solar radio flux appears to be related to how they track thl 
num ei o sunspot groups and the number of individual sunspots observable on the Sun’s disk with 
sunspot number being more heavily influenced by the number of groups and with the . 0 7 cm set’ radio 

th I othdr g h m T ^ y ^ thC nUmbCr of observed suns pots [308], Geomagnetic activity on 

the other hand does not generally peak until after sunspot maximum (during the decline of the sunspot 

the cl'*] /r Ug i 1 , 1S at e evated leveIs froni about 2 to 3 years into the cycle until about 9 to 10 years into 
ve-ir y !fte CyC e " WaS UnUSUal m that ’ while a minimum in the geomagnetic indices occurred about I 
y a ter sunspot minimum, the “true” minimum occurred in 1980, or near sunspot maximum- 
geomagnetic maximum occurred in 1982, being the largest ever recorded in terms of the aa index)' 
Because the geomagnetic indices may be related to the solar wind speed |48, 67-68], solar wind speed is 

speed' sill r C gK Tt' °" UVe , rage ' attCr SUnSPOt maximUm <P robab| y associated with the corotating, high- 
.peed, solar wind streams that are spawned from coronal holes [147, 325)). 

10 7 "Z ha r ng n eSta K iShed ,hat the SUnSp0t cyde Can serve as P roxy lor the solar cycle (and the 
compare ZlZfFlZ ^ “ d S ° lar irradiance >’ with *e caveats as noted, one can 

enorh n E a 1 Van ° US markers of the suns P ot c y de (e.g., average sunspot number 

epochs of minimum and maximum, rise and fall, etc.) to search for possibly significant statistical asso- 
ciations. From Figure I, one finds that more EI Nino events (ignoring intensity) had onsets when the 
annua, sunspot number was below average (54.3) than when sunspot number was above zveLTinZ- 
t.cular, .4 out o 38 El Nino events had onsets when sunspot number was below 54.3, and of the ^4 
events associated with below average annual sunspot number, 17 were of moderate intensity (Table 5). 

TABLE 5. NUMBER OF EL NINO EVENTS IN RELATION TO THE MEAN ANNUAL 
SUNSPOT NUMBER (= 54.3), BASED ON THE MODERN ERA 
SUNSPOT OBSERVATIONS (1848-1987) 


El Nino (Intensity) 

Below (<54.3) 

Above (5*54.3) 

Total 

W/M 

3 

3 

6 

M 

7 

3 

10 

M + 

7 

1 

8 

(Subtotal) 

07) 

(7) 

(24) 

S 

4 

4 

8 

s + 

0 

2 

2 

vs 

3 

1 

4 

(Subtotal) 

(7) 

(7) 

04) 

Total 

24 

14 

38 
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While true, these results should not really be surprising, because for HO years ol the 140-year modern era 
sunspot record the annual sunspot number was below 54.3. Theretore, by chance, one inters that __ 
events should have had onsets when the annual sunspot number was below 54.3 (ignoring intensity), ol 
which 14 should have been of moderate intensity. From the binomial formula (Lapin 1 155], p. 163) and 
using a trial success probability of 80/ 1 40 ( = 0.571), one easily calculates the probability of obtaining 
the observed results, or those more suggestive of a departure from independence (chance), to be P[>24) 
= 27.9 percent and P[ 2= 1 7 1 = 12.3 percent, respectively, and by hypothesis testing, because P|>24] 
and PIS 5 17] are both >10 percent, one infers that the association between the number ol onsets ol El 
Nino and below average annual sunspot number and between the number of onsets ot moderate El Nino 
and below average annual sunspot number, in each case, is statistically unimportant (i.e., the dilletence 
between 24 and 22 events and between 17 and 14 events can, in both cases, be attributed entirely to 
chance). 

The above analysis was based on a comparison of numbers ot El Nino events and the average 
annual sunspot number. Instead, if one uses for comparison the “median" annual sunspot number ( = 
45.3), implying a trial success probability ot 0.5, then one linds that 23 El Nino events had their onsets 
when annual sunspot number was below the median value annual sunspot number, ot which 17 were ol 
moderate intensity. From the binomial formula, one now calculates P[= ; 23] = 12.8 percent and P[> 17] 
= 3.2 percent, the first of nearly marginal statistical significance and the second statistically signilicant 
at >95 percent level of confidence. Comparing the numbers of onsets of El Nino (by strength) against the 
median value of annual sunspot number (i.e. , the lasher’s test for a 2 x 2 table; Everitt |65|. p. 15), one 
finds that a marginally significant association appears to exist between the strength of El Nino and sun- 
spot number, with moderate El Nino favoring the time when annual sunspot number is below its median 
value and with strong El Nino tavonng the time when annual sunspot number is above its median value 
(P = 8.8 percent). Table 6 and Figure 6 summarize these results. 

In Table 7, the number of El Nino onsets is determined relative to whether the sunspot cycle is 
rising (minimum to maximum) or tailing (maximum to minimum). From Table 7, one finds that 27 out ol 
38 El Nino events had onsets during the falling (or declining) portion of the sunspot cycle. By chance, 
one would have expected 24 events ( = 38 x 7/1 I ; a sunspot cycle, on average, lasts 1 1 years, taking 4 
years to rise from minimum to maximum and 7 years to tall from maximum to subsequent cycle 
minimum). Applying the binomial formula, one computes P(5 : 27| = 21.8 percent, interring that the 
association between the number ot El Nino onsets and the phase ol the solar cycle (in particular, the 
falling portion of the sunspot cycle) is statistically unimportant (i.e., the difference between 27 and 24 
events can be attributed entirely to chance). 

In contrast to that reported above, a possibly significant statistical association is found between 
the number of strong El Nino and the falling portion of the sunspot cycle. From Table 7, one sees that 1 2 
out of 14 strong El Nino events had onsets during the tailing portion ol the sunspot cycle. By chance, one 
expects only 9 events. Applying the binomial formula, one computes P|>12] = 6.9 percent, a mar- 
ginally significant result (i.e., at about the 93 percent level ol confidence, one infers that the difference 
between 12 and 9 events cannot be attributed to chance). 

Table 8 compares the occurrences of El Nino onsets to the occurrences of sunspot minimum and 
maximum (actually, a window of ± 1 year bounding each sunspot minimum and maximum). One finds 
that during the 3-year window of sunspot minimum, 1 1 El Nino ol moderate intensity had their on- 
sets. as did two strong El Nino (one of which was classified “very strong"); during the 3-year window 
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TABLE 6. NUMBER OF EL NINO EVENTS FOR ABOVE/BELOW MEDIAN 
ANNUAL SUNSPOT NUMBER (= 45.3) INTERVALS (1848-1987) 


Interval 

Years 

Duration 

Above/Below 

Moderate 

Strong 

Total 

1 

1848-52 

5 

A 

1 


1 

2 

1 853-57 

5 

B 

2 


2 

3 

1858-62 

5 

A 

1 


1 

4 

1863 

1 

B 




5 

1864 

1 

A 


1 

1 

6 

1 865-68 

4 

B 

2 


2 

7 

1 869-73 

5 

A 


1 

i 

8 

1874-80 

7 

B 

2 

1 

3 

9 

1881-85 

5 

A 


1 

1 

10 

1886-91 

6 

B 

1 

1 

2 

1 1 

1892-95 

4 

A 




12 

1896-1904 

9 

B 

2 

1 

3 

13 

1905-08 

4 

A 

2 


2 

14 

1909-14 

6 

B 

1 

1 

2 

15 

1915-19 

5 

A 

1 

1 

2 

16 

1920-25 

6 

B 

1 

1 

2 

17 

1926-29 

4 

A 




18 

1930-35 

6 

B 

1 

1 

2 

19 

1936-41 

6 

A 

1 

1 

2 

20 

1942-45 

4 

B 

1 


1 

21 

1946-51 

6 

A 

1 


1 

22 

1952-55 

4 

B 

1 


1 

23 

1956-61 

6 

A 


1 

1 

24 

1962-65 

4 

B 

1 


1 

25 

1966-72 

7 

A 


1 

1 

26 

1973-77 

5 

B 

1 


1 

27 

1978-84 

7 

A 


1 

1 

28 

1985-87 

3 

B 


1 

1 

Totals 

Above 

70 

14 

7 

8 

15 


Below 

70 

14 

17 

6 

23 


Combined 

140 

28 

24 

14 

38 


24 


EL NINO STRENGTH 

STRONG MODERATE 


DC 

111 

CD 


O 

CL 

C/) 

z 

Z) 

(/) 


CO 

in 

tT 

AJ 


CO 

in 

V 


8 

7 



6 

17 


;z> P[2l7] = 8.8% 


„ , A , V i contineencv table showing the association between annual sunspot number and 

S Irnber^El Nmo events, dividing the events according to strength (modern, e/strong 

and median value (45.3,. The probability of obtaining the — 
suggestive of a departure from independence, is P|^17| - p 

Fisher’s exact test for 2x2 tables. 


° c^POT OBSERVATIONS (1848-1987) 


El Nino (Intensity) 

Rise 

Fall 

Total 

W/M 

0 

6 

6 

M 

5 

5 

10 

M + 

4 

4 

8 

(Subtotal) 

( 9) 

(15) 

(24) 

S 

0 

8 

8 

S + 

0 

2 

2 

VS 

2 

2 

4 

(Subtotal) 

( 2) 

(12) 

(14) 

Total 

1 1 

27 

38 
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TABLE 


8. NUMBER OF EL NINO EVENTS IN RELATION TO THE EPOCHS OP 
SUNSPOT MINIMUM AND MAXIMUM FOR CYCLES 10 TO 2? 


- - — * — ■ i V/ 1 


Minimum ( ± 

1 year) 

Maximum ( ± 1 year) 


Cycle 

Moderate 

Strong 

Moderate Strong 

Total 

10 

1 


1 

2 

1 1 

2 


1 

3 

12 


KVS) 

1 

2 

13 

* 



* 

14 

i 


1 

2 

15 

i 


1 1 

3 

16 

i 



1 

17 


1 


1 

18 

i 



1 

19 

i 


I 

2 

20 

i 



] 

21 

i 



] 

22 

(Subtotal) 

i 

0 1*) 

( 2 ) 

* * 
( 3 *) ( 4 *) 

I* 

Total 

13 * 


7* 

20 * 


NOTE: A moderate event is in progress during Cycle 13 at minimum (having be«un 2 vears 
pnor to minimum); Cycle 22 maximum has not yet occurred Y 
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near sunspot minimum is P1^I3| - ^ a poss ibly significant statistical 

FI Nino onsets near sunspot minimum is I |^l 1 1 ^.u P erce • J >qa 

a s n S ion l found for moderate El Nino events to have their onsets near sunspot mm, mum (at *% 

' ' . , | confidence) Usin" only El Nino events with onsets in the windows ol sunspot minimum 

•*>• one finds a tendency for moderate^ El 

sunspot minimum and strong El Nino to associate with sunspot maximum (P - 7.8 percent, 

A cycle-by-cycle listing of El Nino onsets (by intensity) is given in Table 9, which also identilies 
the maximum amplitude (sire) of the cycle in terms of the annual sunspot ^ 

into two equally sized groups (removing Cycles 9 and 22 because they are only partially complete) 
c'ciITh) to I 5 and 1 6 to 21 one finds that the average number of El Nino events (ignoring intensity) has 
decreased Cycles It) to 15 averaged 3.50 El Nino onsets per sunspot cycle (having a standard deviant n 
of 0.55). while Cycles 16 to 21 averaged 2.33 F.1 Nino onsets per sunspot cycle 

lion of I 03) Hypothesis testing, using the t statistic for independent samples suggests that the di 
ence'in !he two means is statistically significant a, *99.5 percent level of coni, deuce. Thus Ite nmttb r 
of El Niiio onsets per sunspot cycle was greater during Cycles 10 to 1 5 (as compared to Cyc . 


EL NINO STRENGTH 


STRONG MODERATE 



CO 

Figure 7. A 2 x 2 contingency table showing the association between sunspot cycle phase and 
number of El Nino events, dividing the events according to strength (moderate/strong) 
and epochs of sunspot minimum and maximum. The probability ol obtaining the 
observed result, or one more suggestive of a departure lrom independence, is P| 3= 1 1 1 — 
7.8 percent. 
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TABLE 9. NUMBER OF EL NINO EVENTS BY SUNSPOT CYCLE NUMBER 

DURING THE INTERVAL 1848-1987 



Note. Modern sunspot observations began about sunspot maximum during Cycle 9- 
consequently, the S+ El Nino ( 1 844-45) near minimum has not been included in this 
summary; also. Cycle 22 (Rmax = ?) is only partially complete. 



when maximum amplitudes tended to be smaller (larger). Application of Spearman’s rank correlation test 
(Lapin 1 1 55 ] , p. 633) suggests that the number of El Nino onsets per sunspot cycle may vary inverse y 
with the size of the sunspot cycle at >95 percent level of confidence, having a Spearman rank correlation 
coefficient r s equal to about -0.58 (implying that about one-third of the variation found in the number of 
El Nino onsets per sunspot cycle may be due to the relative size differences of sunspot cycles; it should be 
noted however that the inferred inverse correlation between number of El Nino onsets per sunspot cycle 
and the size of the sunspot cycle may be due simply to the “bimodality” of the sunspot cycle, a preferen- 
tial division of the sunspot cycle into two separate groups according to period length: short-period cycles 
with a period <1 i years and long-period cycles with a period >1 1 years, with short-period cycles, on 
average, being about 1.3 years shorter in duration than long-period cycles; ct. Rabin, Wilson, and Moore 
|224|t Wilson 1 300-302. 304|). The implication of this result (i.e., the inverse correlation) is that, 
because of the long-term secular increase in sunspot number (Wilson [303]), more El Nino events may 
have occurred during the Maunder minimum, a protracted period (1645 to 1715) of low sunspot number 
(Eddy [60-62]), than has been experienced of late, and that fewer El Nino events will occur in future 
sunspot cycles (modulated by the 80- to 100-year period or Gleissberg cycle; Kopecky [145]) and pro- 
vided that the secular trend remains upward. From Quinn et al. [221], one finds some support lor an 
inferred increase in El Nino activity during the Maunder minimum in that nine strong events occurred 
between 1652 and 1715. this number of strong events being about 50 percent higher than found for 
today’s cycles. (Moderate events were tabulated by Quinn et al. only for the post- 1800 time span.) One 
should also note that if the height of the tropopause really varies with the sunspot cycle, being lower neai 
sunspot minimum than near sunspot maximum, then volcanic eruptions, perhaps even less violent ones 
of VEI =£4, may more easily penetrate into the stratosphere, thereby initiating climatic change. This 
would seem to be especially so during prolonged periods of low sunspot activity like the Maunder 
minimum than during times of enhanced solar activity. Perhaps this may explain the greater frequency of 
strong El Nino during the Maunder minimum and of El Nino in general during the period of sunspot 
minimum that is seen. 


D. Looking Ahead to the “Next” Anticipated El Nino 

From a statistical point of view, El Nino appears, on average, about once every 3 to 4 years, 
having a range of 1 to 8 years. Assuming that the elapsed time between successive onsets of El Nino is 
normally distributed with a mean of 4 years and a standard deviation of 2 years, one anticipates an El 
Nino to occur usually within 2 to 6 years following the last known onset of El Nino, certainly within 8 
years. Thus, looking at the 1982-83 El Nino as an example, one finds that an El Nino of moderate to 
stronger strength was expected about 1978 to 1982 from the 1976 last known occurrence of El Nino. 
From this perspective, one views the 1982-83 El Nino as an “expected” event. Similarly, having 
observed the onset of an El Nino in 1 982, one expected the next El Nino to occur about 1 984 to 1 988 and, 
indeed, a moderate El Nino had its onset in 1987 (actually, in late 1986; Kousky 1146]). Marking the last 
known Ei Nino onset as having occurred in 1987, one anticipates that another El Nino of moderate to 
stronger strength should make its appearance about 1989 to 1993, certainly by 1995. 

A comparison of El Nino onsets in relation to major volcanic activity showed that about 40 per- 
cent of the El Nino events were preceded (within 1 year) by a major volcanic eruption in the tropics and 
that nearly 70 percent were preceded (within 1 year) by a major volcanic eruption somewhere in the 
world. If the amount of time between eruption (cause) and onset of El Nino (effect) is as long as 2 years, 
the inferred association appears even stronger: about 60 percent for major tropical eruptions and about 80 
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percent lor any major eruption. Thus, a violent eruption of a volcano somewhere in the world often is 
found to herald the occurrence of a moderate to stronger El Nino (an exception may have been the 
1972-73 El Nino, which followed the eruption of Femandina by about 4 years). Such may have been the 
case lor both the 1982-83 and 1987 El Nino events, in that both were preceded (within 2 years) by major 
volcanic activity. It follows, then, that if in the future, a major eruption takes place somewhere in the 
world, a moderate to stronger El Nino might be expected to follow within 2 years. Because major 
volcanic eiuptions occur cju ite often, about once per 1 .6 years ( = 152 years/94 eruptions), in general, or 
about once per 2.1 years (= 152 years/73 grouped eruptions) for “grouped” (same year) eruptions, 
presuming that a major eruption occurs in 1989, one may anticipate the occurrence of the next moderate 
to strong El Nino during the period 1989 to 1991 . (The range of elapsed times between major eruptions, 
based on Table 2, is 0 to 10 years.) 


A comparison of El Nino onsets to various aspects of the solar/geomagnetic cycle also suggests a 
means whereby one might estimate the likelihood ol an impending El Nino. For example, nearly two- 
thirds of the El Nino events had their onsets when sunspot number was <54.3 (the average annual sun- 
spot number, based on the 1848 to 1987 observed values). Further, of the 24 that occurred when annual 
sunspot number was below 54.3, 17 were of moderate strength (or about 70 percent). Likewise, 27 
(about 70 percent) El Nino events occurred after sunspot maximum during the decline of the cycle (when 
geomagnetic activity tends to be greatest) and 12 out of the 14 reported (about 85 percent) strong El Nino 
occurred post sunspot maximum. Regarding the epochs of sunspot minimum and maximum (± 1 year 
windows), one linds that about one-third of the El Nino events occurred near sunspot minimum (13 
events, ol which 1 1 were of moderate intensity) and that about 20 percent occurred near sunspot 
maximum (7 events, of which 4 were ol strong intensity). The 1982-83 El Nino occurred about 6 years 
into Cycle 21 , occurring during the same year as the peak of the geomagnetic activity (which, incidently, 
was the highest peak ever recorded, based on the aa geomagnetic index; the record high monthly value 
occurred in February 1982), and the 1987 El Nino occurred 1 year after sunspot minimum for Cycle 22. 
The year 1988 (2 years into the cycle) saw annual sunspot number surpass the average annual sunspot 
number value ol 54.3, and it should remain above the average value until about 1993 (7 years into the 
cycle); peak sunspot activity is expected to occur about 1989 to 1991 . Thus, the likelihood of another El 
Nino is, at present, increasing and will be especially high after sunspot maximum and also after a return 
to below average sunspot number. From Table 9, one recalls that, for Cycles 10 to 21, every sunspot 
cycle had at least one El Nino, with 3 to 4 events per cycle during Cycles 10 to 17 and only 1 to 2 events 
per cycle since Cycle 18 (which began in 1944). Cycles of late have been among the greatest on record 
(Cycles 19, 21 , and 18 are the three largest sunspot cycles of the modern era, given here in descending 
order). It is now evident that Cycle 22 will also be a larger than average sized cycle (Wilson (304-307]). 
Therefore, one expects fewer than three El Nino events to occur during its course, of which one has 
already occurred (the 1987 event). It follows that only 1 to 2 events may remain to be seen during the 
current sunspot cycle, projected to end about 1996 to 1998. (Also from Table 9, one sees that every cycle 
during Cycles 10 to 2 1 , except Cycle 1 8, had at least one strong El Nino; hence, because fewer than three 
El Nino events are expected for Cycle 22 and because one has already occurred, being of moderate 
strength, the next anticipated event may be of strong intensity.) 

Figure 8 provides a visual summary for much of that just discussed, plotted in the form of 
histograms and plotted relative to the start of the sunspot cycle. On the bottom is a histogram of the 
onsets of the 39 El Nino events appearing in Table 1 , with filled areas referring to strong El Nino events. 
Above the bottom panel is a histogram of the 94 occurrences of major volcanic eruptions (from Table 2), 
where filled areas refer to major eruptions in the tropics (the asterisks serve to note that the mystery 
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events of 1982 and 1983 have been included as tropical eruptions). Above this is a histogram of the 73 
‘grouped occurrences of major volcanic eruptions, where a “grouped” eruption refers to same year 
multiple eruptions. Across the top are major milestones during the sunspot/geomagnetic cycles, with Fm, 
Rm refemn 8 t0 minimum values of the radio, geomagnetic, and sunspot indices, respectively, 
and FM, GM, and RM relerring to their respective maximum values. A casual comparison of El Nino 
onsets and volcanic eruptions (the bottom two panels) suggests striking similarity, especially if one incor- 
porates a lag of 1 year between volcanic eruption (cause) and El Nino onset (effect); in particular, peaks 
of El Nino activity usually follow by 1 year peaks of volcanic activity. Hence, volcanic activity may be 
the dominant causative factor in climatic change, although the sunspot/geomagnetic cycles seem' to play a 
partial role as well. As a future activity, it would be interesting to compare monthly or seasonal values of 
solar/geomagnetic data with respect to El Nino onsets. 


V. SUMMARY 


The major findings of this investigation can be summarized as follows; Onsets of moderate to 
stronger El Nino have occurred at the rate of about one every three to four years, on average, having a 
range of one to eight years. Moderate events have occurred more frequently than strong events and there 
has never been a single occurrence of two strong events in succession. Elapsed times between successive 
onsets may be described using the normal distribution, although a positive skew is apparent with greater 
variation since about 1917 than before. Onsets of El Nino usually occur during southern hemispheric 
summer/fall (December-May) and, as yet, none has occurred during southern winter (June-August). 
Major, tropical volcanic activity (assuming that they inject large quantities of aerosols into the 
stratosphere that persist up to 3 years) can only account for about 70 percent of the El Nino events, 
although major volcanic activity, in general (ignoring latitude), can “explain” about 85 percent of the El 
Nino events. For some unknown reason, moderate El Nino events appear to preferentially occur in the 
vicinity of sunspot minimum and when annual sunspot number is low, especially when it is below 45.3 
(the median value), while strong El Nino events appear to preferentially occur during the declining por- 
tion of the sunspot cycle. A possibly significant inverse correlation is inferred to exist between the 
number of El Nino events per sunspot cycle and the size of the sunspot cycle, suggesting that during the 
Maunder minimum (1645 to 1715) more El Nino events took place than take place today and that during 
the next several sunspot cycles (assuming a continued increase in sunspot number with time, modulated 
by the Gleissberg cycle) fewer events per sunspot cycle might be expected. Finally, one infers that 10.7- 
cm solar radio the aa and Ap geomagnetic indices, and solar irradiance are all directly correlated 
with sunspot number, although the geomagnetic cycle usually lags the sunspot cycle by I year and the 
maxima for the other indices sometimes differ markedly, with geomagnetic maximum occurring about 6 
years into the cycle (or about 2 years past sunspot maximum during the declining portion of the sunspot 
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